Diluted one-dimensional spin glasses with power law decaying interactions 
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We introduce a diluted version of the one dimensional spin-glass model with interactions decaying 
in probability as an inverse power of the distance. In this model varying the power corresponds 
to change the dimension in short-range models. The spin-glass phase is studied in and out of the 
range of validity of the mean-field approximation in order to discriminate between difi'erent theories. 
Since each variable interacts only with a finite number of others the cost for simulating the model is 
drastically reduced with respect to the fully connected version and larger sizes can be studied. We 
find both static and dynamic evidence in favor of the so-called replica symmetry breaking theory. 
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Mean-field spin glass models are known to have a 
rather complex low- temperature phase [ij , which has not 
been clearly observed so far in numerical simulations of 
finite-dimensional models with short range interactions. 
Theories alternative to the mean-field one have been pro- 
posed 0], but short-range systems with quenched disor- 
der are very tough to study analytically Numerical 
simulations have been, thus, extensively employed, devel- 
oping more and more refined algorithms over the years, 
though with no conclusive indication on the nature of the 
spin-glass phase in finite dimension. 

Long-range models are such that their lower critical di- 
mension is lower than the one of the corresponding short 
range model. In particular, one can have a phase transi- 
tion even in one dimensional systems, provided the range 
of interaction is large enough. One dimensional spin glass 
models with power-law decaying interactions actually al- 
low to explore both long- and short-range regimes by 
changing the power [3, [1^0, [3| . These models would be 
perfect candidates for comparing the spin glass phase in 
and out of the range of validity of the mean field approxi- 
mation. Unfortunately, since each variable interacts with 
all the others, numerical simulations are very computer 
demanding and it is hard to get a clear numerical evi- 
dence supporting a specific spin glass theory [^, 0| ■ 

We, therefore, introduce a diluted version of the model, 
where the mean coordination number is fixed (see also 
Ref. H). In diluting, the run time grows as the size N 
of the system, rather than proportionally to N"^ . This 
is a fundamental issue because finite volume effects are 
strong in these models: previous studies were restricted 
to A'^ < 512, while we can now thermalize systems up to 
N = 16384, thus keeping these effects under control. 

We are interested in analyzing the difference among 
the predictions on the spin glass phase of the droplet the- 
ory the trivial-non-trivial (TNT) scenario Q and the 
replica symmetry breaking (RSB) theory [ij. Studying 
the thermodynamics, we focus on site and link overlaps, 
providing strong evidence that both fluctuate in the infi- 
nite volume limit. From the dynamic behavior we learn 



that the four-point correlation function goes to zero at 
large distances when extrapolated at inflnite times. In 
this framework we are able to identify a characteristic 
length-scale £(T; t) for such a decay. 

The model investigated is a one dimensional chain of 
N — L Ising spins (ct^ = ±1) whose Hamiltonian reads 



n = 



(1) 



The quenched random couplings are independent and 
identically distributed random variables taking a non 
zero value with a probability decaying with the distance 
between spins di and Uj, r^j = |i — j| mod {L/2), as 



(2) 



Non-zero couplings take value ±1 with equal probability. 
We choose an average coordination number z — 6 and 
periodic boundary conditions. 

The universality class depends on the value of the ex- 
ponent p, and it turns out to be equal to the one of the 
fully connected version of the model, where bonds are 
Gaussian distributed with zero mean and a variance de- 
pending on the distance as Jf^ oc 0, 0, H, 0] • The 
overline denotes the average over quenched disorder. 

As p varies this model is known to display different 
statistical mechanics behaviors. For the diluted case they 
are reported in Tab. HI In the limit p — *■ the model is a 
spin- glass on a Bethe lattice [lol . 11|, at variance with the 



fully connected version where this limit is ill-defined for 



P< 1 

1< p< 4/3 
4/3 < p < 2 
p = 2 
p > 2 



Bethe lattice like 

2"^* order transition, mean-field (MF) 
2"^^ order transition, infrared divergence (IRD) 
Kosterlitz-Thouless or T = phase transition 
no phase transition 



TABLE L From infinite range to short range behavior of the 
spin-glass model defined in Eas. (|ll2|l . 
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any p < 1. If the decay is gentle enough (p < 4/3), the 
mean- field (MF) approximation is exact. As it becomes 
steeper (p > 4/3), the MF approximation breaks down 
because of infrared divergences (IRD). For p = 2 one 
does not expect a finite temperature phase transition, 
though power law correlations might still be there [2^. 
This special case deserves further investigation. 

The p = 4/3 case corresponds to the upper critical 
dimension of short-range spin-glasses in absence of an 
external magnetic field {D — 6), whereas p ~ 2 plays 
the role of the lower critical dimension. An approximate 
relationship between p and the dimension D of short- 
range models can be identified as follows. In long-range 
models, the free theory in the replica space reads 



n = 



dk 
2^ 



E 

a=ib 



Qab{k) 



(3) 



where a and b are replica indices and Qab{k) is the 
Fourier transform of the distance-dependent overlap ma- 
trix element Qab{i'ij)- Comparing the critical scaling 
(mo cx |r - Tel = 0) of Eq.jll) with that of the free 
theory for short-range spin glass models in D dimensions 
{H ^ J d^k k^TiQ^) the following equation turns out 
to hold close to the upper critical dimension 



1 + 2/D. 



(4) 



We study the equilibrium properties of the diluted 
long-range model both in {p = 5/4) and out (p — 
3/2, 5/3) of the MF regime for sizes up to i = 16384 = 
2^^. We simulate two replicas cr|^'^'' using the parallel 
tempering algorithm [l^ with 20 temperatures and we 
measure site and link overlaps, respectively defined as 
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-T 

N ^ 
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as well as the correlation length [l3| : 



1 



2sin(7r/L) 



Xsg 



X(27r/L) 



(5) 



(6) 



where Xsg = L{q'^) is the spin glass susceptibihty (((• • • )) 
denotes the thermal average and (• • • ) denotes the aver- 
age over the disorder) and x{k) is the Fourier transform 
of the four-point correlation function (x(0) = Xsg)- Av- 
erages over the disorder are taken on 0(3 x 10^) samples 
in the smallest lattices and over 0{2 x 10'') samples in the 
larger ones. In order to compute critical temperatures, 
critical exponents and the Finite Size Scaling (FSS) cor- 
rections we have used the quotient method [iJI • We have 
computed the exponent u using the scaling of the tem- 
perature derivative of ^l/L and rj from the scaling of Xsg- 
As a typical case, we show in Fig.[T]the temperature and 
size dependence of Xsf and ^l. In the quotient method 
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FIG. 1; p = 3/2, IRD regime. Plot of L^'^Xsg vs. T. Inset: 
^l/L vs. T. Sizes are L = 2", with /t = 8, 10, 12, 14. 





P 


Tc 


1/u 


V 


V (th.) 




MF 


5/4 


2.191(5) 


0.28(2) 


1.751(8) 


1.75 


0.40(2) 


IRD 


3/2 


1.758(4) 


0.25(3) 


1.502(8) 


1.5 


0.60(6) 


IRD 


5/3 


1.36(1) 


0.19(3) 


1.32(1) 


1.33 


0.8(1) 



TABLE II: Estimates of critical temperature and exponents. 



the estimates of the critical exponent still depend on the 
lattice size: the extrapolation to infinite volume provides 
both their asymptotic values and the to exponent of the 
leading FSS correction, 0{L^'^). 

The results are summarized in table |lTl The rj expo- 
nent coincides with the theoretical prediction 77 = 3 — p 
(77 is not renormalized in the IRD regime 0,3). Due 
to strong finite size effects this check failed in previous 
works [7|. The ly exponent is consistent with the theo- 
retical prediction, v — l/(p — 1), in the MF case. In the 
IRD regime, thermodynamic fluctuations dominate and 
a renormalization is necessary: at present only one loop 
calculations are available 0, 0] , but the estimate of u is 
too rough to compare with numerical data. 

In the spin glass phase (T < Tc), site and link overlap 
distributions, P{q) and Pi{qi), can be used as hallmarks 
to discriminate among different theories for finite dimen- 
sional spin glasses. Indeed, three cases are contemplated 
in the literature. 

1. Droplet theory: one state; both distributions delta- 
shaped. 

2. TNT scenario: many states {q fluctuates), but 
droplet-like excitations {qi fluctuations vanish for 
large sizes); non-trivial P{q) and trivial Pi{qi). 

3. RSB theory: many states with space-filling excita- 
tions; both distributions broad. 

Distributions P{q) and Pi{qi) for T ~ 0.4 Tc are plot- 
ted in Figs. [2] and [3] in a model case where MF is exact 
(p = 5/4) and in a IRD case (p = 3/2), respectively. In 
both cases, we see two peaks in the Pi{qi) for large sizes 
(insets). Such a result would have been impossible to 
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p = 5/4 (MF) 
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FIG. 2: P{q) and Pi(qi) at p ^ 5/4 (MF). T = 0.9 ~ OATc^ 
with L = 2" and «: = 7, . . . , 13. Inset: Pi{q) for L = 2". 
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FIG. 3: P{q) and P,{qi) at p = 3/2 (IRD). T = 0.7 ~ OATc^ 
with L = 2" and «: = 7, . . . , 12. Inset: Pi{q) for L = 2^^. 



FIG. 4: Distributions of (line) and gaux (empty squares) for 
p = 3/2, L = 2^^ and T = 0.7. Inset: A and B vs. L-°-'^'^ ob- 
tained by measuring the Kullback-Leibler divergence between 
the two distributions for L = 2", k = 6, 8, 10, 12. 



(KL) divergence [18l| between the distribution of qi and 
that of Qaux- In Fig. m for L — 2^^, we compare the opti- 
mal distributions, which should coincide if Eq. ([7]) held. 
Eq. ([7]) provides a strong evidence for a non-trivial link 
overlap distribution as long as the B parameter converges 
to a non zero value for large size, as one can verify in the 
inset of Fig. [?] where we show such an extrapolation for 
A and B plotted vs. an inverse power of L (C and the 
optimal KL divergence go to zero, as expected). 

As a complementary approach we look at the off- 
equilibrium four-point correlation function 



(8) 



observe out of MF with sizes smaller than L = 2^^ [2J]. 

Both distributions seem to be broad, but their ther- 
modynamic limits must be carefully controlled. While 
it is easy to convince oneself that P{q) is not a bimodal 
distribution as L — *■ oo (e.g., P(0) is practically size in- 
dependent), the limit of Pi{qi) is more difficult to extract 
from finite size data, since its variance converges to a 
small value. We provide, thus, an alternative method of 
analysis testing the hypothesis that both q and qi are 
equivalent measures of the distance among states [l5| . 
The simpler relation between the two overlaps is 



gaux = A + Bq^+ C^/T 



q^z , 



(7) 



where z is a Gaussian random variable with zero mean 
and unitary variance mimicking finite size effects, and 
A, B and C are fitting parameters. Such a relation 
is satisfied in the Sherrington-Kirkpatrick model, with 
A = C = 0, and it is a good approximation for the short- 
range spin glass in _D = 3 [3,[l3|- For each value of L, at 
p = 3/2 and T = 0.7, we compute the best fitting param- 
eters by minimizing the symmetrized Kullback-Leibler 



For very large distances the fastest decay expected goes 
like a;"'', because of long-range interactions. For interme- 
diate distances, up to an effective crossover length t{t), 
we observe a slower decay with < a < p, which 
is incompatible with the onset of a plateau at g|.^ in the 
large times limit. This suggests to use the function 
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\1) 



^1/5 



(9) 



to interpolate Cq{x,t) data at a fixed time t. The fits 
are very good and their quality can be appreciated in 
Fig. Ofor an IRD {p = 3/2) system of size L = 2'^\ The 
crossover length £ plays a role similar to the correlation 
(or coherence) length in short range spin glasses 19|, 2^, 
21[. We allow the fitting parameters A, a and d to depend 
on time. Nonetheless we observe (see inset of Fig. [6]) that 
they become stationary for large times: this is a strong 
evidence that Eq. ([9]) is a significant and robust behavior. 

The growth of £{t) with time at different tempera- 
tures below Tc is plotted in Fig. [51 The length £{t) 
reaches very large values (> 10^) with respect to pre- 
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FIG. 5: Cq{x,t) for p = 3/2, L = 2" at T = 0.7 ~ OATc and 
different times t = 2*,2",2i**. The curves are fits to Eq.®. 
Inset: data collapse by plotting C{x;t) x" vs. x/£{t), for all 
times between 2* and 2^^. 



onset of a plateau at any large x. These observations 
are consistent with the clustering properties of the RSB 
theory. The bond diluteness of the model under inves- 
tigation strongly reduces simulation times and allows to 
thermalize systems of sizes large enough to clearly dis- 
cern the double peak structure of Pi{qi)- Both droplet 
and TNT proposal appear, in conclusion, not consistent 
with a FSS analysis over large sizes and with the behav- 
ior of the four-point correlation function and the related 
coherence length. This model is well suited to study the 
spin glass transition in a magnetic field and work is in 
progress in this direction. 

We thank Silvio Franz for useful discussions and sug- 
gestions. This work has been partially supported by 
MEC, contracts FIS2006-08533-C03 and FIS2007-60977. 
Part of simulations were performed in the BIFI cluster. 
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FIG. 6: Crossover length l{t) vs. ^/Tlogt for p = 3/2, L = 2'" 
and r between 0.7 ~ 0.47^ and 1.4 ~ 0.79rc. Inset: parame- 
ters A, a and 5 vs. time at T = 0.7. 



viously studied spin glass models ESUH- In this re- 
gion t{t) is very well fitted by the phenomenological law 
a(T) exp(6(T)-yr log t) , with a and b not very dependent 
on the temperature; this seems reasonable since in acti- 
vated processes the typical scaling variable is Tlog(i). 

We also tried to fit the previous £{t) data with a gener- 
alized droplet scaling [21| t{£) = A(T)€^<= exp(T(T)^'^), 
where the power-law factor dominates near the transi- 
tion (lim^^Tc TC"(T) = 0) and the exponential term gov- 
erns the low temperature regime. The critical exponents 
Zc and -0 are predicted not to depend on T. The data 
shown in Fig. [S] are not compatible with this scaling law 
for any temperature-dependent A{T) and T(T). 

In conclusion, we have introduced a model which is 
easy to simulate and allow to probe the spin glass phase 
beyond mean-field. In this regime, from the analysis of 
thermodynamics, we observe that both the site and the 
link overlap parameter fluctuate for large sizes. In the 
large times limit the out-of-equilibrium four-point func- 
tion Cq{x, t) tends to a well defined function that displays 
a power-law decay to zero and is incompatible with the 
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